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We demonstrate the mixing of rotational states in the ground electronic state using microwave 
radiation to enhance optical cycling in the molecule yttrium (II) monoxide (YO). This mixing 
technique is used in conjunction with a frequency modulated and chirped continuous wave laser to 
slow longitudinally a cryogenic buffer gas beam of YO. We generate a measurable flux of YO below 
10 m/s, directly loadable into a three-dimensional magneto-optical trap. This technique opens the 
door for laser cooling of molecules with more complex structure. 

PACS numbers: 37.10.Mn, 37.10.Pq, 37.10.Vz 


The field of cold and ultracold molecules is undergo¬ 
ing rapid growth with the development of new techniques 
and applications [T] . Cold polar molecules are continuing 
to improve tests of fundamental symmetries [Hg. The 
strong dipolar interaction present in polar molecules also 
yields a rich set of applications in novel quantum mat¬ 
ter [IHS] and cold chemistry [101412) . 

In many of these applications it is desirable to have a 
slow source of molecules to either enhance the interaction 
time or for loading into a trap. The magneto-association 
and coherent state transfer of ultracold alkali atoms [T3] 
can produce polar molecules near quantum degeneracy. 
Supersonic beams produce molecules with forward veloc¬ 
ity > 300 m/s [H], whereas cryogenic buffer gas beams 
have velocities of 50-200 m/s, depending on source con¬ 
figuration [TS] and with the lower velocity beams directly 
loadable into a trap m- Many techniques to further slow 
these beams have been developed. Stark [171 HI] and 
Zeeman [HI [20] deceleration have been used to load con¬ 
servative traps, which can be used for evaporation mi- 
centrifugal slowing [55] and opto-electric cooling was 
demonstrated with CH3F molecules [ 53 ]. The rovibra- 
tional branching that prevents optical cycling transitions 
in molecules was addressed theoretically [241 El] and op¬ 
tical Doppler cooling was subsequently demonstrated in 
the SrF system [H]- Magneto-optical trapping was pro¬ 
posed for TiO |25j and demonstrated in two dimensions 
for YO [ 57 ] and in three dimensions for SrF [531 12H| . 

While it is possible to load directly a buffer-gas cooled 
atomic beam into a three-dimensional (3D) magneto¬ 
optic trap (MOT) [30], for molecules with much lower 
photon scattering rates longitudinal optical cooling is 
necessary for loading into a 3D MOT. Typically 10^ pho¬ 
tons must be scattered to slow molecules to within the 
capture velocities of a normal MOT. During the slowing 
process, the molecular transition used for optical cooling 
experiences a Doppler frequency shift of tens or hundreds 
of resonance linewidths, resulting in inefficient photon 


scattering for a monochromatic unchirped laser beam. To 
date, molecules that have been laser cooled all have the 
magnetic sub-levels in the ground electronic state con¬ 
tinually remixed by the multi-leveled optical cycling pro¬ 
cess. Hence, the widely used atomic Zeeman slower m, 
which relies on the Zeeman shift of a single magnetic sub- 
level to compensate for the Doppler shift during deceler¬ 
ation is not applicable. To maintain a sufficiently large 
optical scattering rate throughout the slowing process, 
two techniques have been employed. Chirping the laser 
frequency was shown to reduce the velocity of a super¬ 
sonic beam of CaF by 30 m/s [35]. For SrF, broadband 
laser radiation was used to slow molecules to be loaded 
directly into a 3D MOT [55]. 

The advent of optical cycling in molecules [SS] [26] is 
critical for the radiation pressure force to directly slow 
and cool molecular beams. However, despite the achieve¬ 
ment of quasi-closed optical cycling in several molecular 
systems, preventing optical pumping into dark states still 
remains difficult for a large class of molecules. While 
vibrational dark states can be addressed with repump 
lasers [53], rotational dark states often have parity selec¬ 
tion rules that prevent direct optical repumping. There 
are several loss mechanisms that break the rotational clo¬ 
sure in the current schemes for molecular optical cycling 
transitions, and they arise from additional decays via 
intermediate states or higher order transition moments. 
Vibrational states typically have lifetimes ~1 s [ 33 ], and 
rovibrational decay is not yet a limiting factor for cur¬ 
rent techniques. Decays via magnetic dipole transitions 
are typically suppressed by a factor of « 5 x 10“® com¬ 
pared to electric dipole transitions (a is the fine structure 
constant), though recent measurements in OH indicate 
that they can sometimes be a factor of 10 faster than 
expected [ 33 ]. Decays through intermediate electronic 
states, present in for example YO [35] and BaF [36], are 
suppressed by the relatively long transition wavelengths. 
For YO, the suppression factor is ^3 x 10“^, however the 
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FIG. 1. (Color online) (a) YO vibronic structure. Dashed ar¬ 
rows indicate decay paths with corresponding Franck-Condon 
factors, q I3Z1. Solid arrows indicate cooling and repump 
laser transitions, (b) Hyperfine structure of the ^4^111/2 

states. The solid arrows indicate the three hyperfine pumping 
components used in each vibrational manifold, (c) Schematic 
of rotational mixing. Dashed arrows indicate the decay paths 
to X^E through the A'^/^ 2,/2 state and black solid arrows 
show the microwave rotational mixing. 


FIG. 2. (Color online) (a) Depiction of the beam apparatus. 
The YO molecular beam interacts with a counter-propagating 
slowing laser beam along x. A 45° probe beam is used for 
Doppler sensitive detection. A PMT is used to collect fluo¬ 
rescence induced by either the slowing or probe beam, (b) 
Two-stage buffer gas cell, (c) Velocity distribution measured 
at the detection region. The solid line is a Gaussian fit with 
center velocity of 69 m/s and full width at half maximum 
(FWHM) of 54 m/s, corresponding to a temperature of 6.7 K. 


decay is sufficient to cause population leakage to rota¬ 
tional dark states that reduce photon scattering. In this 
Letter, we report microwave mixing of rotational states 
to close these additional loss channels. This enables us to 
slow longitudinally a beam of YO emitting from a two- 
stage buffer gas cell to molecular velocities < 10 m/s, 
which can be loaded into a 3D MOT. 

The optical cycling transition in YO proceeds between 
the X^S —>■ A^ni /2 states as described in Ref. [27]. In 
brief, highly diagonal Franck-Condon factors limit vibra¬ 
tional branching m such that a single repump laser lim¬ 
its branching loss to 3.6 x 10““^ (Fig. [^a)). Rotational 
branching is suppressed by a combination of parity and 
angular momentum selection rules [251138| when driving 
the X^E,X" = 1 level, where N is the rotation quan¬ 
tum number. Figure [^b) shows the hyperfine structure 
of the X^E states, where G is the coupled nuclear and 
electronic spin, and F is the total angular momentum. 
For the excited A^ni /2 and states, J is the to¬ 

tal angular momentum excluding nuclear spin. For all 
three electronic states, p = ± represents the parity. 

Figure[^a) shows the additional loss channel where the 
A^XVi/ 2 , J' = 1/2 state decays to the A'^A 3 / 2 , J' = 3/2 
state. Decays to states of higher J are forbidden due to 
angular momentum selection rules. This state has a ra¬ 
diative lifetime of ^1 /iS [35] and will rapidly decay back 
to the X^E state. Since this is a three-photon process 
back to the ground state, parity selection rules will only 
allow decays to even rotational states (Fig. [^c)). In the 
X^E state, while J is not a good quantum number, each 
state of N” > 0 can be expressed as a superposition of 


states with J = X ± 1/2 (for N” = 0, J = 1/2). The 
AJ = 0,±1 angular momentum selection rule prevents 
decay to N" > 4 and thus leads to optical pumping into 
the N" = 0, 2 states. The rotational constant, B, for the 
X^E state is 11.634 GHz. Microwave radiation tuned to 
2B and AB can mix the N” = 0, 2 levels respectively with 
the N” = 1 level. This mixing effectively removes the ro¬ 
tational dark states; however it also lowers the maximum 
optical scattering rate as the number of states in X^E in¬ 
volved in optical cycling increases by a factor of two. 

We employ this rotational state microwave mixing 
scheme for longitudinal slowing of a YO beam generated 
from a two-stage cryogenic buffer gas beam source m- 
We laser ablate (with a 10 Hz repetition rate) sintered 
Y 2 O 3 pellets in a copper cell filled with helium buffer gas 
at 3.5 K (Fig. [^a)). The hot YO thermalizes rotation- 
ally and translationally to the cold helium buffer gas and 
is hydrodynamically extracted through a 3 mm aperture 
into the second stage cell. By virtue of the 4.5 mm gap 
between the two stages and the larger exit aperture di¬ 
ameter of 9 mm, the helium density is ^ 10 times smaller 
in the second stage than that of the first stage and the 
molecules experience only a few additional collisions be¬ 
fore exiting the second stage cell. This slows the 120 m/s 
beam from the first stage m to 70 m/s (Fig. gc)). The 
YO beam then travels 89 cm to the detection region, 
during which it may be illuminated by a 6 mm diam¬ 
eter, counter-propagating slowing beam. The 614 nm 
v” = 0 laser has 70 mW of power and the 648 nm v” = 1 
laser has 80 mW (Fig. [^a)). The slowing laser has three 
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FIG. 3. (Color online) (a) Optical cycling enhancement due 
to rotational mixing. The fluorescence is induced by the high 
power optical slowing beam and v” = 1 repump: (i) No mix¬ 
ing; (ii) IV" = 0 N” = 1 mixing only; (iii) N” = 2^ N" = 
1 mixing only; (iv) simultaneous N” = 2 N" = 1 and 
N" — 0 N" = 1 mixing. The inset shows the fractional en¬ 
hancement of optical scattering vs. microwave detuning from 
the rotational resonance, (b) Comparison of branching losses 
between v" — 1 and (i) Microwaves and v” = 0,1 

are on throughout; (ii) v” = 1 repump is shut off from 14-19 
ms after ablation; (iii) Microwave mixing is switched off from 
14-19 ms. The inset shows curves (ii) and (iii) divided by 
curve (i). Solid green lines are exponential fits. 


frequency components to address the hyperfine manifold 
(Fig.[^b)). To provide the rotational mixing, we use two 
pairs of microwave horns strategically located along the 
beam path, each pair producing mW of either the 
23 GHz or 46 GHz microwaves. 

For velocity sensitive detection, we measure the fluo¬ 
rescence induced by a low intensity 4 mW/cm^, 614 nm 
beam aligned 45° to the molecular beam. This probe 
laser has a saturation parameter of ^0.5 and two fre¬ 
quency components separated by 780 MHz (Fig. [^b)). 
This eliminates mechanical and optical pumping effects 
from the detection beam. The detection laser can be 
scanned over a range of ~100 MHz at a rate of 1 kHz, 
corresponding to a velocity detection range of ~100 m/s. 

We study the effects of rotational mixing by observing 
the fluorescence induced by the slowing beam in Fig.|^a). 


We apply the microwaves in four binary combinations. 
For case i, no microwaves are applied. In contrast, when 
microwave mixing is applied to both N" = 0, 2, (curve 
iv), there is a factor of four increase in fluorescence. Ap¬ 
plying only the N" = 0 ^ N" = I oi N" = 2 ^ N" = I 
mixing (Fig. [^b) ii, iii respectively), we see a factor of 
two increase in fluorescence compared to the no mixing 
case. We estimate the vibrational lifetime of the v” = 1 
state from molecular parameters calculated in Ref. [31] to 
be ^600 ms, which is much longer than the optical pump¬ 
ing time scale. The branching ratios for magnetic dipole 
transitions from the A^ni /2 state is 2/3(l/3) for the 
N" = 0(2) levels in the state, while electric dipole 
transitions from the A'^A 3/2 state have equal branching 
ratios to the N" = 0, 2 rotational levels, consistent with 
curves ii and iii. Hence, decays through the intermediate 
A'^A 3/2 state are identified to be the dominant process. 

To verify that the mixing microwaves are sufficiently 
strong to address all hyperfine states without the need of 
additional frequency sidebands, we simultaneously vary 
the frequency of the N" = 1 o N" = 0(2) mixing 
microwave by 5(25) and measure the corresponding en¬ 
hancement in fluorescence (Fig. [^a) inset). This en¬ 
hancement follows a Lorentzian lineshape with a FWHM 
of 24 MHz, which is larger than the hyperfine variations 
of <5 MHz among the three relevant rotational mani¬ 
folds [101 m] . The mixing for higher vibrational states is 
suppressed because the rotational constant B for v" = 1 
is 50 MHz smaller than that for v” — 0. This results in 
a 100 (200) MHz detuning for the N” = 0 (2) transition 
in v” = 1. Hence, we do not suffer any further reduction 
in the maximum optical cycling rate. 

Figure l^b) shows a measurement of the branching loss 
from the A^Iii /2 state to the A'^A 3/2 state. Again we 
measure the fluorescence induced by the high power slow¬ 
ing beam. For case i, the slowing and repumping light as 
well as microwaves are turned on at all times. For case ii, 
the v” = 1 repump laser is rapidly shut off in ~10 /iS from 
14 ms to 19 ms after ablation. This leads to an abrupt de¬ 
crease in the fluorescence as the YO molecules are rapidly 
pumped into the v" = 1 state and hence turn dark. To 
extract this pumping rate, we divide curve (ii) by curve 
(i) and fit an exponential with a constant background 
ratio of 0.25, resulting in a 1/e time constant of 620 /iS. 
Combining this with the calculated Franck-Condon fac¬ 
tor (7oi = 3.8 X 10’ -3 [37], yields a photon scattering rate 
of 4.3 X 10® s“^. This is smaller than the maximum pos¬ 
sible scattering rate by a factor of ^5 [27] and may be 
due to inefficient depletion due to highly Doppler detuned 
molecules, insufficient power in the wings of the slowing 
beam or an overestimate of the goi Franck-Condon factor. 
Curve iii shows the fluorescence when only the rotational 
mixing microwaves are turned off between 14 — 19 ms. 
The 1/e decay due to pumping into N" = 0 and 2 dark 
states has a 7.4 ms time constant (inset curve iii). By 
combining the ratio of the two loss rates and goi, we esti- 
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FIG. 4. (Color online) (a) Timing sequence used to slow 
molecules to below 10 m/s. (b) Example of spectrally 
broadened slowing light corresponding to a center velocity of 
50 m/s. The black lines shows the sidebands generated from 
an EOM and rapid polarization switching. The red curve 
shows molecular deceleration with an optical saturation pa¬ 
rameter of 3. (c-f) Comparisons between no slowing in red 
and with slowing in blue. Solid lines are guides to the eye. 
(c) Doppler spectra for a 10 ms long slowing sequence, (d) 
Time of flight for 30 m/s molecules, (e) Doppler spectra for 
a 24 ms long slowing sequence as depicted in (b). (f) Time of 
flight for 10 m/s molecules. 


mate the branching ratio from A^ni /2 to to be 

~3 X 10“"*^. This branching ratio agrees with the value 
derived from the transition dipole moments calculated in 
Ref. |39j and is similar to the combined branching ratio 
to v" > 2. At 19 ms after ablation, the vibration repump 
(curve ii) and mixing microwaves (curve iii) are turned 
back on and the fluorescence immediately recovers. 

We estimate the capture velocity of a 3D MOT for YO 
with a laser beam diameter of 1 cm to be 10 m/s. No 
detectable population of 10 m/s molecules exists in the 
buffer gas beam source. To slow some fraction of the 
molecules to below this capture velocity, we apply longi¬ 
tudinal radiation pressure slowing with protocols shown 
in Fig. I^a). To increase the fraction of molecules ad¬ 
dressed by the slowing beam, we spectrally broaden the 
light [32] via a resonant electro-optic modulator (EOM) 
with a modulation frequency of 10 MHz and a phase 
modulation index of 3.3. This spectrally broadens each 
laser frequency to ±30 MHz, corresponding to a velocity 
spread of ±20 m/s. Polarization switching used for desta¬ 
bilizing Zeeman dark states [27] further adds 5 MHz side¬ 


bands and thus we obtain frequency components spaced 
by the natural linewidth of the main cooling transition, 
ensuring that we have continuous deceleration across this 
velocity range (Fig. Qb)). 

As the range of the slowing force is smaller than the 
velocity spread of the YO beam itself, we additionally 
implement frequency chirping on the slowing light. In 
Fig. I^c) and (d), the slowing sequence is designed to 
enhance molecules at 30 m/s. The molecules propagate 
freely for the first 12 ms after ablation, then the slow¬ 
ing beam is turned on from 12-22 ms. During this time, 
the center of the broadband radiation is linearly swept 
from being resonant with 64 m/s to 44 m/s molecules. 
In panel (c), we see a molecular population enhancement 
between 20 — 40 m/s (negligible for below 20 m/s), and 
depletion for velocities greater than 40 m/s. In panel 
(d), we maintain the probe laser’s frequency to be res¬ 
onant with the 30 m/s molecule and we observe a clear 
enhancement in the number of molecules in this velocity 
class compared to the original beam. Due to off-resonant 
excitations, some faster molecules (35 m/s) contribute to 
the detected signal at 27 ms for the unslowed case. 

To produce even slower molecules, we increase the time 
during which the slowing laser is turned on and perform 
a frequency chirp to address lower velocity molecules. 
The slowing light is turned on 2 ms after ablation and 
is left on for the next 24 ms. During the hrst 20 ms, 
the center frequency of the broadband slowing laser is 
linearly increased from being resonant with 64 m/s to 
34 m/s molecules (Fig. |^a)). Figure |^e) shows the ve¬ 
locity distribution after the slowing sequence compared 
to the original velocity distribution. Slowing enhances 
the molecular population below 20 m/s in the detection 
region and depletes molecules with velocities higher than 
20 m/s. Again, we tune the Doppler-sensitive probe laser 
to be resonant with the 10 m/s molecules and we see a 
clear enhancement in the number of molecules in this ve¬ 
locity class (Fig. I^f)). Unslowed molecules leaving the 
cell with 10 m/s forward velocity will arrive 89 ms after 
the ablation, well after the detection time. Hence, we at¬ 
tribute the non-zero 10 m/s signal in the unslowed case 
to off-resonant excitations from fast molecules. 

We have demonstrated enhanced optical cycling in YO 
via microwave mixing of rotational states in the ground 
electronic state. This enabled us to slow longitudinally 
the YO beam via radiation pressure, and hence produce 
molecules that can be loaded directly into a 3D MOT. 
To further close the vibrational branching to < 10“®, 
an additional v' = I ^ v" = 2 repump laser at 649 
nm can be used |271I37|. By turning off the microwave 
mixing for N" = 10 N" = 0 at an appropriate time, 
we can accumulate trapped cold YO in the ro-vibrational 
ground state. The A'^A 3/2 state lifetime is ~1 ^s, which 
is >10 times longer than that of A^IVi/ 2 - This opens 
up the possibility of a narrow line MOT smij with a 
Doppler temperature of 10 /rK as compared to the 116 /rK 
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obtainable from the d' = 0 to transition. 
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